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Abstract 

The real Ginibre ensemble refers to the family ofnxn matrices in which each 
entry is an independent Gaussian random variable of mean zero and variance one. 
Our main result is that the appropriately scaled spectral radius converges in law to a 
Gumbel distribution as n — > oo. This fact has been known to hold in the complex and 
quaternion analogues of the ensemble for some time, with simpler proofs. Along the 
way we establish a new form for the limit law of the largest real eigenvalue. 

1 Introduction 

In [T5] J. Ginibre introduced the basic non-Hermitian ensembles of random matrix theory. 
These are n x n matrices M comprised of independent (and standardized) real, complex or 
quaternion Gaussian entries, and are clear analogues of the Gaussian Orthogonal, Unitary, 
and Symplectic Ensembles (G{0/U/S}E). 

The results of [15 include explicit formulas for the joint density of eigenvalues of M, 
in both the complex and quaternion cases. The real-entried case posed serious technical 
hurdles, due largely to the fact that the real line itself receives positive mass in this setting, 
and the determination of the joint eigenvalue density remained open until |10j |23j . These 
papers found conditional densities for the real Ginibre ensemble eigenvalues, that is, formulas 
for the joint law given a predetermined number of real eigenvalues. Even with these in 
hand, the expressions proved sufficiently complicated that the finite dimensional correlation 
functions — the basic tool(s) required to obtain limit theorems for local eigenvalue statistics 
— were only determined in the last few years. In [3] Borodin-Sinclair rigorously established 
that the eigenvalues of the real Ginibre ensemble form a Pfaffian point process: there are 
2x2 skew matrix kernels for the real/real, complex/complex, and real/complex correlations 
from which the general fc-point correlation is built as a 2k x 2k Pfaffian. These formulas 
along with the connected skew orthogonal polynomials discovered by Forrrester-Nagao [T5] 
allowed Borodin-Sinclair to derive the scaling limits for the kernels at both the (real and 
complex) bulk and the (real and complex) edge [5]. (Concurrently, Sommers |30) reported 
the scaling limits of the kernels in the bulk) . 

Here we are after a scaling limit for the spectral radius in the real Ginibre ensemble. 
Among other motivations, this may be viewed as one possible refinement of the circular law. 
The latter refers to the fact that the normalized counting measure of the scaled eigenvalues 
converges (weakly almost surely) to the uniform measure on the unit disk. This result has 
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a rather long history, starting with the work of Girko |16j which was made rigorous by Bai 
PQ, and culminating in the universality (in terms of entry distributions) theorems of |18] 
and then [33]. On a local scale, considerable progress has been made on the universality of 
the n f oo bulk correlations (even all the way up to the edge) for both real and complex 
entries, see [51 El [35] 

Our main result reads: 

Theorem 1.1. Denote by R n the spectral radius of the real Ginibre ensemble and set j n = 
log(n/(27r(log n) 2 ). Then, as n — > oo, 




where G is the Gumbel law with distribution function Fc(t) = e ? e . 

One can certainly adjust the scaling so that the limiting distribution function takes the 
more standard form of e~ e . It is written this way for comparison: at the same scaling the 
limiting spectral radius in the complex Ginibre ensemble is also Gumbel, with distribution 
function e~ e . A similar result holds in the quaternion case. The universality of the limiting 
Gumbel law for spectral radius in any setting (real, complex, or quaternion) has not been 
addressed. 

The analog of Theorem |1.1| for complex Ginibre is a triviality. The eigenvalues of the 
complex Ginibre ensemble form the canonical radially symmetric determinantal process on 
the complex plane, and it is the case that the moduli of the points of such a process are 
independent. Stated in this generality this fact can be found in [HI Ch. 7], but had been 
observed earlier by Kostlan [2U] specifically for complex Ginibre. While a Gumbel law 
is the only possible scaling limit for the extremal point, the precise scalings were worked 
out in [26] where the author was unaware of Kostlan's result but rediscovered and used a 
consequence thereof. In the quaternion case there is nothing like this determinantal trick. 
Still the expectations of certain eigenvalue class functions factor nicely, and this turns out 
to be enough, see again [26 . 

For real Ginibre, there appear to be no shortcuts toward pinning down the fluctuations 
of the spectral radius. Instead we return to the standard random matrix theory machinery 
of tracing through the limiting Fredholm determinant /Pfaffian formulas for the related gap 
probabilities (which are available once the correlation functions are known). The proof of 
Theorem |1.1| follows from determining the real and complex gaps separately: 

Theorem 1.2. Let z±, . . . , z n be the eigenvalues of M , then 




max \zk\ < „ 

k:z k eC/R V 4 V^Jn, 

for any t € K as n — > oo, where again j n — log(n/(27r(logn) 2 )). 
Theorem 1.3. Introduce the integral operator T with kernel 

T(x,y) = ~ e- (x+u)2 e~ {y+u)2 du. (1.1) 

I" JO 

Let x be the indicator of {t, oo). Then 

max z k < y/n + t ] = y/dct(I - T X ) T t , (1.2) 

k:z k £M J 
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where T t is built as follows. Set g(x) — -^j^e x I 2 , G{x) — g{y)dy and denote by R(-, •) 
the kernel of the resolvent operator (I — T X )~ l . Then 

T t = (1 - at) t)te) +l^-h)J (J - T X )- 1 g{x)dx, 

for a t = J7° G(x)(I - T X )- 1 g{x)dx and b t = (I - T X )- l G{t). 

In particular Theorem |1.1| just recasts the result for the complex points, the largest real 
eigenvalue simply lives on a smaller scale. (Obviously the largest negative real eigenvalue 
exhibits the same limit law.) From a technical standpoint this means that we never have to 
deal explicitly with the mixed real/complex correlations. 

A convergence result along the lines of Theorem 1.3 was presented in [13]. There the 
right hand side of ( 1.2 ) was left in terms of the Fredholm Pfaffian of a 2 x 2 matrix operator. 
Here, besides rigorously establishing the appropriate norm convergence, the factorization 
in terms of scalar operators coincides with the initial form of f3 = 1 Tracy- Widom (TW\) 
distribution function as originally found for GOE. Indeed, the resulting structure of the 
above limit law is precisely the same as that form of TW\'. T replaces the Airy kernel and 
the Gaussian density g(-) plays the role of the Airy function Ai(-). Unsurprisingly, our 
derivation of (1.2 1 follows [33] quite closely. 

One point of interest is that T, like the Airy or related Bessel operator, is product 
Hilbcrt-Schmidt, and so trace class, on any positive half-line. (Additional properties of T 
are needed to show that the factors which make up T t are sensible). Unlike the Airy or 
Bessel case however T does not posses a Christoffel-Darboux form, and so is not integrable 
in the sense of [9] , or at least not in this simple way. This presents at least one roadblock in 
obtaining a "closed" expression for ( 1.2 ), say something in terms of a single special function 
like the Painleve formulation of the Tracy- Widom laws. Even a characterizing PDE or 
system of ODEs has eluded us. The full large deviations of (1.2) are also open. The right 
tail, as t — ¥ oo, is easily seen to have a Gaussian shape (this again was pointed out in |13|). 
The determination of the left tail lies deeper and will be pursued in a later paper. This 
seems worthwhile given the separate interest in the limiting largest real point due to its 
applications in the stability analysis of certain biological systems [24) . 

The next section recalls what is needed of the real Ginibre correlation functions and gap 
probabilities. The limit law for the largest complex and real points are derived separately 
in Sections 3 and 4. Section 5 reports on some numerical simulations we have performed 
and discusses additional open questions. 

Acknowledgements. Many thanks to Torsten Ehrhardt and Harold Widom for useful dis- 
cussions. The first author was supported in part by NSF grant DMS-0645756 as well as a 
grant from the Simons Foundation. The second author was supported in part by NSF grant 
DMS-0801243. 



2 Determinants 

The results of [5] lead to (Fredholm) determinantal formulas for the relevant gap probabilities 
in the n x n real Ginibre ensemble M. We state things in the particular form that we need. 

Proposition 2.1. The probability Pc,n(t) of there being no complex eigenvalues of M of 
modulus greater than t > is given by 

P c , n (t) 2 =det(I-K nX ) (2.1) 
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in which K n is a 2 x 2 operator defined on L 2 (C+) © L 2 (C+) ; C+ = {z : Im(z) > 0}, cut 
down by the indicator function X = X{\z\>t}- I n the standard notation, 



K„ 



S n DS n 



with the various operators defined most easily through their kernels: 

ie~^ z ~^ 2 _ - 
S n (z,w) = -= (w-z)4>(z)4>{w)e zw e n ^ 2 (zw), 



/2n 

DS n {z,w) = —iS n (z 1 w) , and IS n (z,w) — iS n (z,w). The shorthand 

— Z K 



) = Verfc^Mz)!), e„(z) = ^-, 



(2.2) 



(2.3) 



fe=0 



is used, where as usual erfc(z) = ^= J z e * dt. 

Note that while [5] considers only even values of n (n = 2M there), the subsequent 
results of 29 shows that the formulas remain unchanged for n odd. This can be loosely 
understood by considering that any "extra" particle must be real. 

For gaps on the real line we have the following. 

Proposition 2.2. The probability of there being no real eigenvalues of M larger than t > 
— oo is given by 

V R .n(t) 2 = det (I-K nX ), 

where again K n is a 2 x 2 operator now defined on L 2 (M.) © L 2 (M.), and x — X[t,oo)- The 
overall form of K n is similar to the complex case: 



K n = 



S n DS 
'IS n 4~ £ S— 



Here e is the operator 



and the basic kernel is 



tf{x) = 2 J s & n (y - x )f(y) d y, 



(2.4) 
(2.5) 



S n (x,y) 



e 2 



v l \ X e 2 

— ;= — e " y t n ^2\ x V) 4 1= 

V2vr 1 ' V^{n-2)\J 

in terms of which DS n = SS^ where 8 acts by differentiation on the first variable, and: 
1. When n is even 



(2.6) 



IS n (x,y) = eS n (x,y); 



2. When n is odd, 

IS„(x,y) = eS n (x,y) + 



1 



eS n -i(x,y) + 



2™/ 2 r(n/2) J 
1 



rx ry 



2ir(n-2)\J Jo 



(w — u)(wu) 



n-2 e -\u 2 -\w 2 dudw ^ 
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In both cases e also acts on the first variable. 

While the first form of the IS n kernel for n odd may be more attractive, the second 
is better for asymptotics. Both structures are also valid for, for instance, for odd n GOE. 
In particular, in the first form the correction term may be expressed as for ir n -i the 
relevant skew-orthogonal polynomial and s n _i its normalizer, see below, and the same holds 
for GOE with appropriate substitutions for 7r„_i and s n -\. 

Note that since the e operator is not trace class, what is meant by the Fredholm deter- 
minant of (2.4) is at first not clear. However, this is a standard technicality in the world of 
f3 = 1 ensembles and how to mollify things — or smooth out the e is well understood (see 
for § VIII]). 

Proof. The n even case is a restatement of [5j Theorem 8] . The n odd case was stated without 
proof in [T^] (based on their derivation of the correlation functions of /3 = 1 ensembles of 
odd order). See also [21] ■ We include the odd n case here since the details of the derivation 
have not before appeared in the literature. We appeal to [3S1 §7] (providing yet another 
derivation of the correlation kernel for j3 = 1 ensembles of odd order) using the particulars 
of Ginibre's real ensemble. 

Recall the weighted skew-orthogonal polynomials for Ginibre's real ensemble [H| given 

by 

n 2j (x) := e-* 2 /V J ; n 2j+1 (x) := e^' 2 ^ 1 - 2jx 2j ~ 1 ), 
with the normalization 

rj := (7r 2j -,7r 2i+1 ) = 2V^(2j)!. 

(We will not explicitly use the skew-inner product, but the interested reader is referred to 
28J). One further set of normalizations (which does not arise for even n) is 



s 2j :-- 



n 2 j(x) dx = y +1 ' 2 T{j + 1/2). 



(In general, similar normalizations are necessary for the odd degree skew-orthogonal poly- 
nomials, but in this case these vanish). In particular, 

Sjl = 2^/ 2 r( J + i/2) = _j_ 

r 3 2>/27r(2j)! ^ +1 jV [ ' ' 

where the last identity uses the Gamma function duplication formula. We will also need 

eir 2j+1 (x) = e- x2 ' 2 x 2 3 and «r n _i(a;) = -T^ 1 sgn(a;) 7 (~, y) ■ 

It suffices to establish (2.6) for n odd, since the DS n and (the modified) IS n terms follow 
immediately from [29]. For S n , [29] p. 31] implies 

c i s e , s 7r n _i(x) ^4 s 2j , . e-K n -i{y) ^4 s 2 j , > n n -i(x) 

S„(x,y) = S n - 1 {x,y)-2— — — ^ -^-e7r 2 j+i(j/) + 2— — - — ^ -^-^j+t {x) 1 

Here, we use 



3=0 3=0 J 



2 E 7^ e ^ + i(y) = e-v 2 ' 2 J2 = e-y 2 /\jMy 2 m (2.8) 

3=0 3 3=0 3 ' 



G 
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and 



to write 



-l(y) _ sgn(y) T (2' V 2 ) 



2 r(f) 



S n (x,y) = S n .i(x,y) + (l - e^e^u 2 ^)) 

- e - 2 / V - 2 sgn(u) 7 (^ 



27r(n- 2)! 



2' 2 



Next, note that 



(2.9) 



2^r(^) Jo 



u n-2 £ -±u 2 d 



more simply, 7 (|, ^) = 2 2 + 1 Jjp u" *e 2 n " ; sothat 



2n{n - 2)! 7 



2ir(n- 2)1 



,71 — 1 „— ^14 



5" dti. 



Since n — 1 is even we can substitute for S n -\ to get 



S n [x,y) 



e 2 



^-y) 2 

-e xv t n -s(xy)+ 
Z7T v27r(n — 3)! Jo 



,71 — 3 _— 4 U 



^(n-2)! 7o 



5" du 



n-2_-iar /•!/ 



n-2 -i« 2 j ^ e 2 



27r(n- 2)! Jo 



u n - l e~^ u du. 
(2.10) 



Integration by parts on the second term yields, 



x n 2 e - 2 x f V u n-3 e -±u 2 du= e fX " e * r (^) r 



x e 2 



V27r(n-3)! J 
The first of these terms takes 



/ 2 7 r(n-2)! V27r(n - 2)! J 



e -i(*-2/) 2 
X/27T 



e -i(x-y) 2 

e xy e„_ 3 (a;u) to -== — e xy e n - 2 (xy), 



'2-k 



and the second cancels the last term in (2.10) to produce the advertised formula. □ 



2.1 From Correlations to Gap Probabilities 

Once again, what [5] and [29] establish are Pfafnan formulas for the fc-point (any combination 
of real and complex) correlation functions. For completeness we briefly review how to go 
from the correlations to the above determinantal gap formulas. 
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The situation is similar to the classical j3 = 1 (GOE) situation, but is complicated by 
the presence of both real and complex eigenvalues. Letting L represent the number of real 
eigenvalues and M the number of complex conjugate eigenvalues, there is a different joint 
eigenvalue density for each pair (L, M) with L + 2M = n. Representing this density as 
SIl,m ■ K L x C M -t [0,oo) (the exact formula for which can be found in [531 US]), the 
normalization constant for the ensemble is given by 



(L,M) 
L + 2M = 



where /iR and /Ltj| are Lebesgue measure on R and M L , and /xc and u.^ are defined analo- 
gously. 

We define the £, m correlation function Re tTn : M. e x C m — » [0, 00) by 
i?,, m (x,z)- — (L-^)!(M-m)!2 M - m 

L>i,M>m 

x / / L . M (xVa,zV/3)d^(a)^- m (/3), 

where, for instance xVa = (x\, . . . ,xe,ai, . . . ,aL-i)- That is, the £, m correlation function 
is a weighted sum of all marginal densities formed by integrating out L — i real variables 
and M — m complex variables from all Ql,m for which this makes sense. 

The main result of [5] demonstrates the existence of three matrix kernels K^' R ,K^' C 
and K^' C (and its transpose K^' R ) such that 



t (x,z) = Pf 



Now consider a C C C which is invariant under complex conjugation, written as the disjoint 
union C = A U B where 4CK and B C C \ R. Conditioning on the the number of real and 
complex conjugate pairs of eigenvalues, we have that 



- C.n 



; P(no eigenvalues in C), 

P(exactly L real eigenvalues and all eigenvalues in C c ) 

L + 2M — n 
1 ^ 1 



Z n ^ L\M\2 M 

{L> ,M) 
L + 2M=n 



x / / (n( 1 -^K))n( 1 -^(^0))^,M(a,/3)d/4(«)^c / (/3), 

with xa and %s the characteristic functions of A and B. Expanding the products in the 
integrand and simplifying leads to: 

= E >J l9 m / / i?^n(x,z)^(x)^(z). 
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In particular, when B = 0, that is when we are interested in the probability that there are 
no eigenvalues in some subset A or K but we place no restrictions on the complex eigenvalues, 



/> i * J q i ■ J A^ 



d/4(x). 



Similarly, if A = 0, 



n/2 

E 

m— 1 

n/2 

E 

m— 1 



(~i) T 

?/;,! 
(-l) r 



Pf 



s,t=l 



(B+)*> 



Pf[^ c (, s ,z 4 )]™^( Z ), 



Where -B + is the component of B which lies in the upper half plane. Each of the last two 
displays defines the Fredholm Pfaffian of the indicated matrix kernel K n , or, in symbols 

® I, and / is the nxn identity matrix. One may then invoke 



Pf(J-iT) where J 



-1 



the relationship between the Fredholm Pfafhan and the Fredholm determinant [5S1 13] : 

det(7 + JK) =Pf(J -Kf. 



3 Complex points 

Returning to \2.2\ and ( |2.3[ ) the correct scaling can be implemented as in: z, w Z n , W„, 
Z n (r,8) = (Vn+ ] f^+-^=y 9 , W n (s, V ) = (Vn +] f^+^=y\ (3.1) 
where j n = log(n/(27r(log n) 2 ) and it is implicit that r, s > —j n - Next, we replace K n (z,w) 

by 

K n (r,0;s,r,) = J |Z "J |Wn| K n (Z n , W n ), (3.2) 

the variable change occurring entry- wise in K n , which acts on 

C t = L 2 (T)®L 2 (T), T= [t,oo) x (0,7r) 3 (r t 6),(s,rj) t > -oo. 

The restriction to T is from now on assumed in the dehnition of the transformed K n . 

The typical procedure (which is followed in the real case) is to identify a limit ker- 
nel/operator K on T for which K n — > K in trace norm, concluding the convergence of 
det(I — K n ). Here instead, though K n is trace class (it is finite rank) it is more convenient 
to cast things in the Hilbert-Schmidt norm, in which K n vanishes in the limit. This prompts 
the introduction of the regularized determinant 



det 2 (J + A) = det((J + A)er A ). 
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In particular, if A is Hilbert-Schmidt with eigenvalues {\k(A)}k>o, there is the evaluation 
det 2 (I + A) = n^LoC 1 + A fe )e" Afc IV.7]) allowing us to write 

det(J - K n ) = dct 2 (/ - K n )e- trRn . 

With a matrix kernel, the trace is just the sum of the traces of the diagonal entries. The 
proof of Theorem 1.2 is then completed via the basic estimate: with 1 1 • 1 1 the Hilbert-Schmidt 
norm, 

|det 2 (Z + A)-det a (/ + B)| < \\A - B\\ exp + + I) 2 ) , (3.3) 

again see [T7l IV.7], along with the next lemma. 
Lemma 3.1. We have that \\S n \\2 — >• as n —¥ oo while 

In addition, \\DS n \\2 = \\IS n \\2 — > 0. norms are with respect to L 2 (T x 7"). 



In particular, the estimate (3.3) gives the desired result upon choosing A — K n and 



B = 0. For a different way to understand Theorem |1.2[ the proof of Lemma [3. 1| will show 
that on bounded sets the kernel S n is well approximated (after an unimportant conjugation) 

by + 

S K (rAs,v) = ^ (1+ e K) 2 ei(e _ r , ) _ 1 , (3.4) 

in which k = n(n) tends to zero as n — ¥ oo. (So, formally, the limit operator has kernel 
2^X{e}( 7 ?) e ~ 3 ' r+s ' ) - Again though, we don't attempt to carry out a proof in this manner.) 
Here the decoupling of the moduli and phases of the points is made explicit, as^is the 
asymptotic independence of neighboring phases (on the scale k). The kernels for DS n and 
IS n will be shown to exhibit a shaper decay. One can also check that det(J— S K ) ~ e _ s e 
from the series definition of the Fredholm determinant. 

As for the proof of Lemma 3.1 we will use the next three estimates on the polynomial 
e«(z). 

Lemma 3.2. (Wimp 0)- Uniformly in t > 0, 

l z n (nt) = l < t <i + -^^crfc(v^A*W) + ° (^)) ' 
where fi(t) = y/t — log t — 1 is taken positive for all t. 

Lemma 3.3. (Bleher, Mallison 0). For small enough 6 > 0, let now fx(z) — \J z — log z — 1 
he uniquely defined as analytic in \z — 1| < 8 with + x) > for < x < S. Then, for 
any M > 1 it holds: with n — > oo, 

e~ nz t n (nz) = — j^- erfc(Vn/i(«)) [l + O 



-tit * 



2y/2n.'(z) V \ n ( z - !) 

for z satisfying ^= < \z — 1| < 6 and |arg(z — 1)| < 
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Lemma 3.4. (Kriecherbauer, Kuijlaars, McLaughlin, Miller [21 ). For any < a < 1/2, set 
U = {\z — 1| < n~ a } and denote by D the unit disk. Then, 



e n -i(nz) = e n z r ' 



i2nrn(l - z) 



1 + 



n|l - zp 



for z€ D c -U 



id all n > 1. 



Both [2] and [21] contain more detailed and complete asymptotics along the lines stated 
in Lemmas |3.3| and |3.4| we record only what is used here. Also, as is easy to check, the 
appraisals of all three lemmas apply without change to e„_2 (rather than say e n ,e n _i) for 
n large enough. Last we should point out that Lemma |3.2| may be arrived at by combining 



Lemmas 3.3 and 3.4 (the asymptotics must be consistent after all). It is convenient though 
to have a single result to quote for the full range of the real argument. 



Proof of Proposition \3 . 1\ This is split into five steps. Throughout, C is a large positive 
constant that may change from one line to the next. 

Step 1 is the trace calculation, integrating 



S n (r,6:r,9) = 



with again Z n = (y/n 



\Zn\ 



2^ 

1 \z n 



: Im(Z„)e 2Im(z " )2 erfc(^Im(Z n )) ( 



Zn-2{\Z n \ 2 ), 



= )e ie over (r,0) G T, recall (3.1|, (3.2). 



Since, for real y > we have that erfc(y) < ^=^6 y while erfc(y) = -^^e v (1 
Q{l/y 2 )) for y -> oo, 



l im -Llm(Z„)e 2Im ( z ") 2 erfc(y2Im(Z„)) = ± 

pointwise on T, with the left hand side being bounded by the right for all large n. 
Next, with 



(3.5) 



tn — —\Z n \ 2 — 1 

n 



\fn An 



Lemma |3.2| implies that 

e -|^l%„_ 2 (|Z„| 2 ) 



1 V(t n )t r . 

sft t n -l 



eTfc(y/rifj,{t n )) 1 + 



(3.6) 



(3.7) 



uniformly on T- And since /u(l + e) = (1 + 0(e)) for < e <C 1, (3.6) and (3.7) produce 

\Z„ 



-\Zn\- 



t n -2{\Z n \ ) = e r (l + o(l)) uniformly for r = 0(^/7^). 



Here we have used that 



(3.8) 
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which in effect dictates the choice of j n . For the tail we have the bounds: with r > and 
n large enough so that j n > 1 while 7n^ 2 /n 1 / 2 < 1/2, 



-\z n \- 



e„- 2 (|^„| 2 ) < ^i 2 cxp(-7V(i„)) 

In 



„ \/n o In 

< C— r 2 exp — 

< cv 2 e - r / 4 . 



n n J \ y/n 



The second line uses the inequality: 



1 



e - log(l + e) > ^(6 - 6 2 )e, for aU e > 5 and < 5 < 1, 



(3.10) 



with the choices £ = t„ — 1 and S = \J^/ n /n. Hence, dominated convergence coupled with 
Q and @ yields 

tr(5 B ) = | S n (r, 9; r, 8)drd9 = V'(l + o(l)), 

as required. 

Step 2 considers S n away (though just barely) from the diagonal. All further nontrivial 
behavior occurs when the argument of z n -2(Z n W n ) is in a small neig hborhood of -Z n W n = 
1 for which we can invoke Lemma 3.3 For given 6 € (0, 7r) consider the set 

■A/e.n = {(»?, r, s) : 77 g (0, tt), |0 - r/l < 1 t < r, s < n 1 / 4 } n T, 



in connection to which it will be useful to make the definition 

' (3-11) 



Similar to before, define z n — A=Z ni w n — -j=W n , and record 



Ffo r + s + r / + s/ 'y/jn) 

n 2^j~n An 



-i(6- n ) 



(3.12) 



Note \z n \ — 1, |u>„| — 1 = 0(e) on and that z n w n satisfy the assumptions of Lemma 

13.31 there: 

;W — < |1 - z n w n \ < Ce, | arg(l - z n w n )\ < tt/2. 
As we still have the estimate erfc(z) = + 0(1/ z 2 )) for arg(z) < 37r/4, it holds that 



-ZnWr. 



z n . 2 {Z n W n ) 



1 



/27m z n w n - 1 



1 + 



1 



since /it' (z) 



z-l 
2zn(z) 



The rational term in (3.131 may be bounded roughly as: 



<CJ— onAf e ,n- 

In 



(3.13) 



(3.14) 
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And combining the exponent in (3.13) with e ^ Zn Wn ^ 2 <fi(Z n )(j)(W n ) , S n has the overall 
exponential factor of n times 



-(zn - w„) 2 -Im(z„) 2 - Im(w„) 2 - \x(z n w n ) 2 

= 1 - \{\Zn\ 2 + \W„\ 2 ) + log(|z„||w„|) 

+ i((Re(z„)Im(z„) + arg(z„)) - (Re(w n )lm(w n ) + arg(u>„)) 



(3.15) 



where we are assuming Im(z„), Im(w„) > 0. The relevant part of (3.15) satisfies 

1 - \(\Zn\ 2 + Kl 2 ) + l0g(K|K|) < - \[{\z n \ - l) 2 + (\w n \ 2 l) 2 ] 



(3.16) 



< 



In r + s 



r 2 + s 2 7; 



1 ITn _ In 

2 V n 4n 

3/2 



2n 



8n 



8j n n n 3 / 2 ' 



for all r, s > and j n /n < 1. Here again (3.10) is used (twice) with e = \z n \ — 1, |s n | — 1 
and 5 = 2V~»F" (For r, s < 0, or in particular just bounded, a Taylor expansion produces 
a better bound — with ^2 replaced by ^rp.) 
Finally there remains the prefactor, 



|Zn||«> n | 



-i(w 



n ^n I 



1 

47rV7n V lHi(2„)Im(u>„) 

1 HT sin(g±2) +3 n i (0 -,) 



(3.17) 



27r V 7n \/sin(0) sin(r/) 
where <5 n is an additive error term that satisfies 8 n = 0(7i _1//4 7n X ' 2 ) = 0(e) for r, s < n 1 / 4 . 



Combining the above (and recalling < 3.9 )), we have the upper bound: 



S n (r,6;s,r)) 



< C 



Sin(^) + e (H-.) 



And since 



•y/sin(0) sin(?7) 

: 1 ,■<>+* ( sin (e+2) +£ ) 



■ Mb. 



(3.18) 



. 2 : ' "/ d V d9 < Ce, 
sm(v) sm(r)) 



(3.19) 

(this for any small e > 0) with a similar bound in a neighborhood of 9 = it, it follows that 



d.O 



< Ce, 



where now recall that e = e n \. 0. 

For the integral over < < e, we go back to the start and bound the erfc(a) appearing 
in each copy of cf>(a) in a simpler way: for a > 0, erfc(a) = -JL J °° e~( t+a ) d£ < e~ a . 
This removes the integrability issues due to the inverse sines in ( 3.18[ ) at the expense of an 
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additional factor of y/n. Importantly though this keeps the over all exponent from (3.15) 
unchanged. In particular, the estimate, 



S n (r,9;s,r}) 



< C\fn I sin 



n 



_ (r+£) 

e I e s . 



is also available on any Ng, n . On the region of current interest sin (j^^j = 0(e) and so 

\S n \ 2 < Cne\ 



JM e 



explaining in part the choice that e should decay a bit faster than n 1 / 4 . 



Remark. If we further restrict r,s = o(y / 7^) the bounds (3.14) and (3.16) can be supplanted 

by 

p i(0-n) 



z n w n - 1 (1 + ^3Ey(fl-.j) - l 



and 



i- 2(W 2 + KI 2 ) +lo g(WKI) 



In _ r + s 
2n 2n 



sin(-p) 
■J sin(0) 8X11(77) 



respectively. Also considering a fixed 9 £ (0,7r) with \6 — 7}\ = one has 

1 + 0(1). 

Therefore, setting f n (r,9) = n(Re(2;„)Im(z„)) + (n + l/2)arg(z„) (recall (3.15)), there 
is the estimate 



1 Ffn 



-(r+s)/2 



2tt V n (1 + ^TS) e i{e-v) _ i 



(l + o(l)), 



in the "bulk" of T, as advertised in (3.4), with n = y^f • 

Step 3 considers again r, s < n 1 / 4 , but keeps e and 77 are away fro m th e diagonal via 
\0 — n\ > e, the latter defined in (3.11 ). With now \z n w n — 1| > ^e, see (3.12 ), we have that 



^(Im( Z „) 2 +Im(«,„) 2 ) e -f( Z „-ffi„) 2 e r lZ „ 1 B„ en _ 2 ( ^^z^ - n ) 



(3.20) 



.1 



1' ,2 1, . ,2 



<C-exp [n(l--\z„\ - -\w n \ + \og\z n \\w n \ 



< c—e~~ [r+s)/8 . 
~ ne 

Lemma |3.4| is responsible for the first inequality, producing the same exponent as in Step 2. 
Line 2 then reuses the estimate (3.16) fro m tha t step in conjunction with (3.9). 

Next, while the considerations behind (3.17) still hold, here the sin(^i 3 ) + e is of little 
use. Instead we are lead to the bound: 

\S n (rA8,ri)\<cJ^~-=== on \9 - V \ > e, r lS <n^ 4 . (3.21) 
\l n e v /sm(6') sm(r?) 

First keeping 9 and 77 away from the origin, let 

O n = {e<9,ri< vr/2, \6 - r?| > e r, s < n 1 ' 4 }, 
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for which we have that 



\SJ'<C- 



^ 2 d9 



sin(0) 



< C- 



1/4 ' 



having substituted (3.11 1. The same bound holds for the integral over the set analogous to 



O n but with 9, rj < n — e. 

To finish, as in step 2 we control the integral over the region where say < 8 < e by 
altering our initial bound on the <j> function(s). Here though this is done in just the variable 
near the singular point. To illustrate, the bound 



\S n {r,e ]S ,ri)\<C 



In n 
n e 



1/4 _(r+«)/8 



= c ln - 



,-(r+a)/8 



is again valid throughout the region described in (3.211, but useful only when 9 is small 
where it produces: 



o n n{e<e} 



S n \'<C ln / \\og(8 + e)\d8<C^ 



7 2 n 



/4' 



Once more, like considerations apply to 9 near tt (and also of course to the situation where 
9 and rj change roles). 

Step 4 dispenses of the case that either r or s is greater than n 1 / 4 . Once again \z n w n — 1| > 



II > Ce, and 



e -n(Im( Z ,,) 2 +Im(«,„) 2 ) e -i( Z „-fi)„) 2 e n 2 „^, en _ 2 ^ ZnWji ) 



1/2 

In 

„3/4' 



< C ln c ~ (r+s)/8 



exactly as in (3.20), now just employing the definition of e. The relevant bound on the 
\S n (r,9;s,r,)\ < Cn^ 4 ln (l + \r\ + \ s \) e -^/ 8 . 



kernel becomes 



Here we have again used the simplified bound </>(a) < e a2 / 2 ; as well as the even rougher 
estimate \z n — w n \ < (C + \r\ + \s\) in the prefactor. In any case it is enough. The square 
integral of the above restricted to {r V s > n 1 / 4 } is dominated by Ce~ n 1 l c . 

Step 5 is to note that everything above applies to DS n (or IS n ) with one notable change. 



The appearance of (w n — z n ) (or (w n — z n )) in (3.17) rather than (w n — z n ), leads to the 
replacement of the factor [sin(^i 3 )+e] with [sin(— ^) +e]. The latter is O(e) on any Af n .e, 
producing an additional decay along the diagonal. This completes the proof. 

□ 



4 Real points 

We run through the calculation over even values of n, returning to the modifications required 
for n odd at the end. 
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4.1 n even 



First the determinant of (2.4) at finite n is reduced to that of a scalar operator. Throughout 
this section any x appearing on its own denotes x = X{t<x<oo}- 



Lemma 4.1. With K n defined in (2.4) we have that 

det(7 - K nX ) = det(J - T nX ) det(J - W n ). 
_ffere T„ is £/ie symmetric part of S n , recall pr 



T n (x,y) 



2tt 



'e„- 2 (a;y), 



and W n is a finite rank operator defined in (4-4)- 



This step in particular mimics Tracy and Widom's treatment of GOE 
Next, introducing the scaling as in 

T n (x,y) = T n (y/n + x,y/n + y), 



(4.1) 



(4.2) 



quite closely. 



(4.3) 



the convergence of the first factor in (4.1), and more, is dealt with by the following (the 
point being that det(7 — T n x) = det(7 — xT n x)- 

Lemma 4.2. For all t > -co, the L 2 operator x^nX converges in trace norm to Tx with 
the kernel forT defined in (1.1). Further, x^nX xT n X an d (I ~ xT n x) 1 (I ~ X^X)^ 1 
in L 1 , L 2 and L°° operator norms. 



The last step deals with the W n operator appearing in the second factor of (4.1). The 
proof of Lemma |4.1| will show that W n is of the form 



W n = a x <8> fa + a 2 ® fa 



(4.4) 



in which 

a 1 = (I - Tnx)" 1 ^ 



a 2 



\ fan, I-XW- T nX Y X K + (I - T nX )- X T n (I - x)) 



fa = X^n 

fa — S t - Soo, 



4>n{x) = Kn J u n ~ 2 e~z u du, and ij) n (x) = K' n x n ~ 1 e~^ x (with certain constants K n ,K' n ). 
The determinant of I— W n is then comprised explicitly of the L 2 -inner products (a*, fa)i<i,j<2, 
and what we need is: 



Lemma 4.3. After scaling as in (4-3), the inner products ({on, (3 
formal limits. 



j)i<i,j<2 



converge to their 



The object identified as T t in the statement of Theorem |1.3| is just the expansion of 
"det(J- Woo)". 

Before the proofs of Lemmas |4.1[ |4.2| and |4.3| we verify, as indicated in the introduction, 
that the kernel T does not have a "Christoffcl-Darboux" structure. Gerard Letac showed 
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us this short argument. The question is whether there exist functions F and G (which can 
assumed to be C 2 ) for which 



e -(*+«)V(y+«) 2 du 



F(x)G{y) - F{y)G(x) 



x-y 



The answer is no. Since the left hand side is of the form e x e y H(x + y), it is enough to 
prove that 

F^GM-F^G^^ 

x-y 

is impossible except for F and G proportional and so H — 0. Making the change of variables 
t — x — y and s — x + y we find that 



31 
dtt 



{F(s + t)G{s - t) - F(s - t)G(s + t)) = 0. 



This implies that 

1 1 — ^ F(s + t)G(s -t)- tF(s + t)G'(s - 1) + tF'(s + t)G(s - t) 

is an even function. Differentiating this function with respect to t and setting t = yields 
F'{s)G(s) = F{s)G'{s) which was the claim. 



Proof of Lemma 4-1 We start with the matrix kernel 

K n = 



tS n + e S n 



(4.5) 



This is exactly analogous to the kernel for GOE given by Tracy and Widom in [3S] and we 
follow the strategy laid out in §11 of that paper. First, since 5e = —I, 



XKnX 



X S 
X 



-6S n x Six 
(-eS n + e)x Six 



Using the famous dct(I — AB) = det(I — BA) trick, we find the Fredholm determinant of 
XKnX equals that of 



-eS n x Six 
(-eS n + e)x Six 



XS 
X 



-eS n xS Six 
-eS n xS + exS Six 



The determinant is further unaffected if we subtract the first row from the second and then 
add the second column to the first, resulting in 

six ~ tSnxS Six 
exS 

The best way to understand this is to note this pair of moves is effected by K n i— > PKn P^ 1 
1 -1 



with P = 



1 



VrW 2 = det 



Thus, 



1 - s lx + tS n xS -six 
-exS I 



det 



1 - s lx + tS n xS - SlxexS 
/ 
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which follows by 'row reducing' the matrix. And since one may check that eS n — S^e, we 
find that for even n, 

P„, R (i) 2 = det(7 - Six + Sn(l - x)exS)- (4.6) 

The above manipulations have been carried out completely formally, with no attention as 
to in which space(s) the operators/determinants reside. The needed technical details may 
be taken (yet again) verbatim from [35] i see §VIII. 
Now for even n recall the kernel 

Sl{x,y) = -L e -^)V^ e », 2 (xy) + l- yn ~ le ~ V J 2 [ X u n - 2 e-^ 2 du 
V2tt V2tt [n-2)\ J 

= T n (x,y) + U n (x,y), 

where we have previously defined T n as the symmetric part of S n , and consider now U n the 
remainder. For later it will be useful to express U n (as an operator) as 



U, 



n ~~ f n 



®i) n = Un j u n - 2 e-^ u2 dv\ <g> (<y"- x e-^ 2 ) , (4.7) 



where K n = ,/ / — -, 1/2 and n' = . " 1/2 - . Keep in mind of course that 6 n — ed>' , 

y v27r(n— A)\ 11 y lityn— 2)\ n 

4>' n (x) = n n x n ~ 2 e~^ x2 . 

Next, introduce the resolvent 

(I - Tnx)- 1 =1 + Rn, or R n = T nX (I - xTnXY\ (4.8) 



with kernel R„(-, •). The determinant (4.6 1 factors as in 

P„jt(t) 2 =det(/-T„ X ) (4.9) 
x det(7 - ((/ - T n x)"Vn) ® (A) + (I - T^SUl - X )e X S), 

having used A(B eg) C)D = (AB) ® (D T C). The second term (in which ex8 appears) is 
simplified by considering the commutator: 

be, <*] = ~( s t ^St-S^® Soo), 

where 5 a is the dirac delta. Since again eS = —I, 

(1 - X )e [ X , S] = (1 - xhxS = -(1 - X)(e* ® S t ~ £oo ® <5oc), 

with now 

e*^) = j s S n ( i ~ a; ) ( so e oo(z) = x)- 

Thus, 

(7 + R n )Sl(I - x)e X S = ~\{I + Rn)Sl(l - X ) ® («t - <U- 

Expanding out the S„ and noting (j4 <£> B)(C ® D) = (B, C)(A (g) D), the second factor in 
(4.9) is the determinant of the identity minus the finite rank operator 

W n = ((I-T nX )- 1 <t>n)® 

+ ^ (W**. 1 - - T n x)-^n + r n x) _1 Tn(l - X)) ® (*t - *oo), 
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to which we apply the well known fact 

det(J - W n ) = det(<y i)j - 

with, as announced above, 

ai = (I - T„xrVn Pi=xip n 



L<ij<2 



«2 2 



* ((V>„, i - x)(/ - r„x) _1 0n + (/ - r n x) _1 T„(i - x)) fa = s t - s 



(4.10) 



Here (aj,/3j) are regular L 2 inner products. (We have reused 5 many times - here it is the 
standard Kronecker delta.) In particular, all components comprising the original det(J — 
K n x) are well-defined. 

We conclude with a few simplifications. First, 

(ai,ft) = (x0„, (J ~ %T„xrVn) (4.11) 
(ai.A) - ((/ - T nX )- x 4>n) (t) - ^n(oo). 

In the second line (7 — T n x)~ 1 <i> n = fin + T n x{I — X^nX) -1< ^n is used. Next, since (T„ + 
T n Rn)x = it follows that 

((/ + J R n )T n (l-x),xV'n) = (l-X,-R„Vn), ((/ + J R„)T n (l-x),5 t ) - (I- X ,Rn(;t)), 
and, with c n = (^„, 1 - x), 

{a 2 ,/3x) = ~(c n (a 1 ,/3 1 )-c n + J (I + R n )^ n (x)dx^j (4.12) 

(a 2 ,j8a) = 2 f c „(ai,/3 2 ) + / R n (x,t)dx 



Combining (4.11) and (4.12), det(I - W n ) equals 

(1- (on, ft)) (l- i / ^(i.tjds^ - J(ai,,9a) / (J" + R n )ip n (x)dx, 



note the c ra factor has dropped out. This last expression has precisely the same structure 
as equation (41) in [35 . □ 



Proof of Lemma \4-ty We employ the trace convergence criteria of Theorem 2.20 of |27j . 
showing that: 

/>oo />oo 

/ T n (x,x)dx —> / T(x,x)dx, (4-13) 

J t Jt 

and 

/>oo />oo pOO fOO 

f(x)T n (x,y)g(y)dxdy-> J J f{x)T{x 1 y)g{y)dxdy, (4.14) 

for all f,g € L 2 ([i, oo)). Both follow from the pointwise convergence of T n to T, and an 
easy domination. 
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Though T n — » T pointwise already demonstrated in [5] , it is useful here to establish some 
local uniformity. Start with the expression 

— c 2 2/) x ~\~ y xy 

T n (x,y) = - — = — e _m7 "e n _ 2 (w7 n ), »fo = Vn(x,y) = 1 + 



Repeating the estimate from Lemma 3.2 (as used in step 1 of Lemma 3.1 ) yields 

T n (x, y) = e ~i iX ~ V)2 ^^eMV^KVn)) (l + O (—)) , 

uniformly, granted that rj n > which holds say for x and y bounded and n large enough. 
Since \fn^(rj n ) — > for uniformly for x and y on compacts, it follows T n (x,y) — > 

e 2 < ^/2^ ) er f c (j^ij m ^ ne same fashion. A change of variables shows this object is equiv- 
alent to T. By the smoothness of the functions involved we also have a constant C so that 
T n (x,y) < CT{x,y) for all x and y bounded. 

For the rest of domination, on diagonal rj n = (1 + x/^/n) 2 is always nonnegative and we 
can continue as above. In particular if x > 1, 

^JheM^iVn)) < r ? n . e-^ < xe~*\ 

Vn - 1 V n {Vn - 1) 

since /i 2 ((l + a) 2 ) > a 2 . This is enough to conclude that (4.131 holds. Off diagonal, let 
x +y > 1) go back to the definition of T n (x, y) and note quite simply that \e~ nrin e n -2(nr) n )\ < 
e~ nr i™+ n \ r >n\ . Hence, when xy > as well we have that T n (x,y) < e~^^ x ~ y ^ which controls 
that range of the integral in (4.141 : J e~^ x ~ v ^ f(y)dy € L 2 for / e L 2 . On the other 
hand, if for instance x > and y < (requiring t < 0), the same obse rvation gives 
T n (x,y) < e~i x e 4 '*^ which suffices for the remaining variable range in (4.14). 

The trace norm convergence certainly implies the L 2 operator norm convergence of x^nX- 
More directly though, any symmetric kernel operator of the form x^X nas L p i— >• L p norm, 
for p = 1, 2, oo bounded as in 



||xM x ||<sup/ \M(x,y)\dx. 
y>t Jt 

(The L 2 i y L 2 bound, less familiar than the L 2 (g> L 2 kernel norm bound, is due to Holmgren 
[2"2l §16.1].) The estimates above will then imply that xT n % -> X T X in the L 1 and L°° 
operator norms too. As for the resolvents, it is enough to check that 



sup / T{x,y)dx < sup I 
y>t Jt y>t V n Jo 



e -(n+yf du _ / e -^du < 1 



for t > — oo. □ 



Proof of Lemma \4--3\ The scaling is to shift t by y/n in all appearances of x (and so R n ). 
This shift then filters into <j) n ,ip n and so on by changing variables in each (a,i,fij). Again 
all scaled functions/operators are decorated with tildes. 

With now g(x) — ~^^ e ~ x ^ 2 > G(x) = g(s)ds one readily finds: pointwise as n — > oo, 

4> n (x) -> #(x), n (x) -> G(x). (4.15) 
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The first convergence also takes place in L 1 n L 2 , while the latter can be considered to hold 
in K + L 2 {s, oo) for any s > — oo upon writing <fi n (x) — (f> n (oo) — J (j>' n . 
Allied considerations will also show that 



/oo />oo 
T n (x,y)dy^> / T(x,y)dy, 
-oo J — OO 



(4.16) 



pointwise and in (L 1 n L 2 )(t, oo). 
Starting with (ai,fa) we find that 

(Si,/3i) = f ^„x, (J - X^nX)" 1 ^™) 

= (0„(oo)x, ^n) ~ ((</>n(oo) - 0„)X, (/ - XTnXy 1 ^)- 

The con verge nce of 4> n {oo) — > 1 can be viewed as holding uniformly, and, by the second 
item in (4.15), 4> n (oo) — <j) n {-) converges in L 2 (t , oo). Furthe r, (I — x^nX) VnG) — > (I — 
xTx)~ 1 9{ ) m L 1 n L 2 by the first item in (4.15) and Lemma 



4.2 



The next three items are treated in much the similarly (to each other). Take J_ R n (x,t)dx, 
and rewrite the scaled version as: 



R n (x,t)dx = / T n x(I - X T nX) 1 T n (x,t)dx+ / T n (x,t)dx 



(4.17) 



The first term is the L 2 (t, oo) inner product of the functions 



T n (x,-)dx and {I - X T nX) 1 T n (-,t), 



(4.18) 



each of which converges in L 2 (t, oo) by (4.16 ). The second term converges to /_ T(x, t)dx 

30 rQ 
-OO Jt 



also by (4.16 ). (As in [35], it is most convenient to see this by writing J* R — (J^° — J^°)R n 



before applying the identity inherent in (4.17)). 

The term f_ (I — T n x)~ 1 ip dx is easier. Now we have that 



(I + T n x) 1 ip n dx= / T„x(I-xTnX) 1 ipn(%)dx+ / ip n (x)dx 



The only real change is the replacement of (/ — %TnX) T n (',t), appearing in (4.17), with 
(I — xT n x)~ 1 ipn(-)- We already have noted that this tends to its formal limit in L 2 . Finally, 
f_ ip n (x)dx is the same as <j) n {t) up to trivial factors and also converges to G{t). 
Returning to (a\, fa), we only need deal with 

/•oo ^ 

T n X{I -X T nX)~ l( t>n(t) = / T n (t, x) (I - X T nX)~ X K (x) dx. 



Again, one can decompose (j> n (x) = </>„(oo) — (</>„(oo) — <t> n ( x )) an d alternatively use the 
L x (t, oo) or L 2 (t, oo) convergence of T n (t, ■) coupled with the L 2 + L°° convergence of (/ — 

XTnXT^nix). □ 
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4.2 n odd 

By Proposition |2.2| when n is odd 



i,„(t) 2 =det(I-K nX ), 



where 



K r , 



-eS n . 



■ e + 



6 ,. 



where 0„ is as in (4.7) and y> = etp n with also -0 as in (4.7) 



Performing the same maneuvers as for the even n kernel, we find that Vs., n (t) 2 (n odd) 
equals the determinant of 

Again we can factor out the (J — T n x), and are left with two extra components as compared 



with (the second factor on the right hand side of) equation (4.9). 

To see that the first extra component, featuring ~ Sn—i> gi ves n0 contribution to the 
determinant in the limit note the following. From before X (T n — 2n-i)x ~* 0: m trace as 
well as the L 1 and L°° operator norms, so this piece may be taken out as a perturbation to 
the (previously factored) (I — T n x). Also, we have that (f> n — n _i — > 0, and ip n — V'n-i — > 
(pointwise and in R + L 2 (s, oo) or L 2 , respectively). The latter two differences will appear 
as factors in some appropriate "a's" and "j9's" in the limiting finite rank determinant, so it 
is enough that all L 2 inner products in which they figure are zero. 

Similarly, writing 



the vanishing of <f> n — tp n pointwise and in M. + i 2 (s,oo) is sufficient to conclude that any 
limiting inner product in which these terms enter in will also be zero. This completes the 
verification. 



5 Numerics and open questions 

Not having a closed form for the limiting distribution of the largest real point prompted us 
to carry out some straightforward simulations of the matrix ensembles, resulting in a few 
notable observations surrounding this object, as well as the finite n behavior of both the 
largest real and (in absolute value) complex points. (From now on we write just "largest 
complex point or eigenvalue" — that we mean in absolute value should be understood.) 

Figure [T] compares the histograms for the largest real and complex points at n — 36,64, 
and 100. Noticeable right away is the heavy left tail in the real point distribution. One 
might have expected that the tail going into the bulk of the spectrum would be lighter than 
that held down simply by the Gaussian weight. On the other hand, recall that there are 
only 0(y/n) real eigenvalues for n j" oo [11], and so rather weak level repulsion along the 
real line. As the right tail of this law can be seen to have Gaussian decay, a reasonable 
conjecture is that the limiting left tail is exponential to leading order. 
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Figure 1: Normalized histograms for the modulus of the largest complex eigenvalue (left) 
and the largest real eigenvalue of 40,000 random n x n matrices for, from lightest to darkest, 
each of n = 36, 64 and 100. 

A closer look at Figure[T]also shows that the empirical distribution of the largest complex 
point appears far more symmetric than its limiting Gumbel shape would suggest. Figure [2] 
focuses in on the n — 100 case, and highlights that at least for moderate n the real point 
distribution is heavier tailed to the right as well. This gets right into some basic questions 
on the speed of convergence for these laws. 




Figure 2: Histograms of the largest complex (dark) and real (light) eigenvalues of 40,000 
random 100 x 100 with an enlarged picture of the right hand tails. 

Experts of RMT anticipate fast convergence of finite n statistics to their limiting distri- 
bution. Here though the Gumbel shape of the spectral radius will not "kick-in" until n is 
considerably large. What dictates the phenomena is that the real point is centered about 
y/ri while the complex point is centered about i/n + yf%x with 7„ = log ( 2 ) ■ Even 
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as written, this is a purely asymptotic statement, j n is not even positive until n ss 165, and 
in any case one sees that n has to be much larger still until the competing real/complex 
distributions separate. 



Figure 3: The proportion of samples for which the largest eigenvalue was real as a function 
of n. The unshaded region was sampled 10,000 times for each n; the shaded region 1,000 
times for each n. 



A first question might then be: what is the chance, at finite n, that the spectral radius 
comes from a real point? Figure [3] graphs this empirical probability, indicating this is a 
slowly decaying function, still just under .4 for n — 100. Quantifying this analytically 
appears challenging. To start, one would need sharp control of the mixed (real/complex) gap 
probability, given by the Fredholm determinant/Pfaffian of a 4 x 4 matrix kernel operator. 

Together with Figure [2j the phenomena appears to be as follows. At moderate n (e.g., 
n = 100), the largest eigenvalue is most likely complex. However, in the situation where the 
largest eigenvalue is real, it is more likely to be larger than if it were complex. With hindsight 
one can see this in the most famous picture attached to the real Ginibre ensemble which we 
repeat an instance of here in Figure [6j The striking feature is the so-called "Saturn effect" , 
based on which alone a person might be forgiven for having conjectured that the largest 
eigenvalue would be real, with probability one, as n f oo. Rather, the Saturn effect is a 
phenomenon which appears from plotting the eigenvalues of many matrices simultaneously. 
Eventually, the complex points overwhelm the 0(y/n) on the real line. 

In summary, one cannot expect the Gumbel law to be a good approximation for the 
spectral radius at small n. And this is due more than just the mixture of the separate 
largest real and complex point laws. Figures [T] and [2] already show that the largest complex 
point distribution itself is not well approximated by its limiting Gumbel. Figure [5] makes 
this more transparent. The issue with ^/7^ not being sensible for smaller n is circumvented 
by appealing to equation (3.9 1, lim„„>— ^ P ~i n / 2 



1, which more or less defines 7„. 
For numerical comparisons then we take 7„ to be the solution 7 of 7 2 e 7 = 

More confounding in Figure [5] is that scaled spectral radius appears better approximated 
by the Gumbel distribution than does the largest complex eigenvalue. This though must be 
purely superficial. It again comes back to the fact that when the largest eigenvalue is real it 
tends to be larger than were it complex, thickening the right tail of this histogram to look 
more Gumbel. As emphasized many times, however, this phenomenon vanishes in the large 
n limit. 
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Figure 4: A plot of the eigenvalues of 1,000 random 100 x 100 matrices. 




Figure 5: A scaled and shifted histogram of the largest complex eigenvalue (left) and the 
spectral radius (right) of 40,000 100 x 100 matrices as compared to the Gumbel density 




In light of all this, a fair question that remains is how to engineer a decent fluctuation 
approximation for the spectral radius at finite n. 

Simpler questions such as determining just how slow the speed of convergence of the 
largest complex point is to its Gumbel limit would also be interesting. Working through 
the proof of Section 3 only produces an 0((log n) -1 ) speed estimate. There is no reason to 
expect this is close to optimal. On the other hand, our last figure compares the spectral 
radius in the n = 100 complex Ginibre ensemble (in which there are no real points with 
probability one) to its corresponding Gumbel limit. The fit is far more satisfying. Studying 

the proof from [25] of this limit theorem gives an 0{ ( lo y) ) speed. 
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Figure 6: A histogram of the spectral radius of 40,000 complex random 100 x 100 matrices 
as compared to the Gumbcl density e~'~ e . 
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